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Одна из ведущих концепций современной теории управления заключается в 
достижении главной конечной цели на каждом этапе функционирования системы, 
которая обеспечивается путем оптимизации объекта управления (ОУ) в реальном 
масштабе времени [1]. Это требует полного использования имеющейся априорной 
информации в виде моделей ОУ и возможно при: 

- оптимальном оценивании (фильтрации) динамических процессов в ОУ; 

- идентификации (оценивании структуры и параметров модели) ОУ; 

- синтезе оптимального управления на каждом этапе функционирования 
системы; 
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- адаптации (настройке оптимального управления при неполной информации). 

Разнообразие и сложность технологических процессов вызывает необхо- 
димость применения универсальных по целям и эффективных за результатами 
принципов управления, что реализуется с помощью оптимального управления. 

Применение метода аналитического конструирования оптимальных регуля- 
торов (АКОР) Летова-Калмана для нелинейных ОУ приводит задачу синтеза к 
поиску решения уравнения Беллмана — нелинейного дифференциального уравнения 
в частных производных относительно вырабатывающей функции, которая 
определяет закон оптимального управления. 

Развитием теории АКОР Летова-Калмана является принцип минимума 
обобщенной работы, разработанный академиком А.А. Красовским, согласно 
которому оптимизация управления осуществляется по функционалу обобщенной 
работы (ФОР) [1, 2]. Преимущество метода АКОР по ФОР заключается в том, что 
его функциональное уравнение, в отличие от уравнения по АКОР Летова-Калмана, 
представляет собой линейное дифференциальное уравнение в частных производных, 
которое имеет принципиально более простые решения. 

Создание систем управления (СУ), как правило, включает следующие этапы: 

- выбор критерия (функционала, цели) управления; 

- разработка математической модели ОУ; 

- синтез законов (алгоритмов) управления; 

- разработка алгоритмов адаптации законов управления по режимам 
функционирования ОУ; 

- реализация полученных законов. 

Для нестационарного технологического процесса такой подход недостаточно 
продуктивен, поскольку структура и параметры алгоритмов управления выбираются 
на этапе проектирования для конкретных (стационарных) условий. Для компенсации 
нестационарности необходимо объединить последние три этапа, что позволит 
создать на базе вычислительных средств управляющую систему, которая 
осуществляет синтез управляющих воздействий и их реализацию в процессе 
функционирования ОУ (объединенный синтез управления). 

Рассмотрим решение задачи оптимального управления по ФОР в процессе 
функционирования непрерывным ОУ с обобщенной дискретной моделью вида: 


ЖЕ+т|= ЕР], иК], Е, С], а[], к}; & =, и-1, 7=0,12..., (1) 

где и - глубина прогноза (для компенсации чистого запаздывания и времени 

на синтез и реализацию управления); Р’ — обобщенная функция (алгоритм) 

преобразования; ХА], и], К], СК], а[к] — соответственно, координаты состояния 

процесса, его управления, возмущения, шума и параметров к текущему такту 

времени К с соответствующими глубинами памяти; /,,К,‚,, — начальные такты 
последовательных этапов (циклов) управления. 

Для реализации оптимального управления необходимо соответствующее 
информационное обеспечение — априорная (адекватные прогнозирующие модели) и 
апостериорная 2[А] (текущие измерения и результаты обработки) информации: 

2] = 4], Ш], А]. (2) 

Для преобразования обобщенной модели (1) к виду с линейно входящим 
управлением выполним расширение пространства управляющих воздействий [1]. 
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Для этого в качестве управления будем использовать изменение значения 
управляющей координаты на текущий такт: 


и [А =и[-иА-П. (3) 
Тогда модель (1) с учетом (3) приобретает расширенный вид с линейно 
входящим управлением и [К] : 


мЕ+п] о Е - 0 Е] 4 
и[К] иА-П| |1 ’ 


или более компактно 

[А+] = А+В-и"[А], (5) 
где Х[А+п]= А+ Пи]; А= {ЕЩЕ-Пу; В= {0,1”. 
Синтез оптимального управления выполним по стохастическому ФОР с 


аддитивными функциями затрат на управление и дискретным временем [1]: 
А! Ка-п-1 Ку 


л=Е, (А, .Р+ >, ОСЬ" [ЕК + 2. И (и ЕТК) + 2. О’ (и[АТЬК)}, (6) 


КЕ, +т 


—п1 


где Е- математическое ожидание; Г,— терминальная ба конечного 


3 


состояния этапа управления (целевая функция); О, (И,- положительно 


определенные функции Ляпунова (функции затрат на качество управления и само 
управление); И,* — положительно определенная функция Ляпунова, которая 


; х29д — заданное значение состояния 


* 
принимает минимальное значение при и = Ио > ь 


процесса; и„ — искомое оптимальное управление, которое доставляет минимум 
функционала. 

В соответствии с принципом разделения, оптимальное в смысле ФОР (6) 
управление стохастическим процессом (1) в условиях некоррелированности целевой 
функции и ошибок измерения ( (7, С) =0) может быть приблизительно получено [1, 


2], как оптимальное управление детерминированным процессом с точным 
измерением вектора состояния 2[А] путем замены его действительного значения на 


оценку по его условному математическому ожиданию 2[А]= Е, {2[]}. Полученные 
таким путем приближенные решения задачи синтеза закона управления 
стохастическим процессом тем точнее, чем выше точность оценивания (меньше 
ошибка |2[А]- 2[%]| ). 


Дискретное уравнение Беллмана с учетом (5) записывается в виде 


И) = ЕСИ, "ПИ, (А+ Ви (П-ОВ, О) 
где = -п-Ь К ыы 


+1 
} 


При предельном условии И, „О =У,О,„]) решение рекуррентного 


уравнения (7) определяется в соответствии с равенством 


д * ОТ. = 
м, Пен В 8 
т (8) 


ор 


и при квадратичной функции затрат на управление равно: 
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ОТ. . (9) 
би +П 

Таким образом, алгоритм синтеза оптимального управления по ФОР с 
прогнозирующей моделью в процессе функционирования системы управления 
включает: 

Г) оценку текущего состояния ОУ в моменты начала очередного цикла 


управления (д, согласно (2); 


*ж . 
И ор [Й Ай 


2) прогнозирование свободного движения ОУ по модели (1) на заданном 
интервале [К, +п, К а —!] оптимизации управления; 


3) вычисление градиента целевой функции Г; +1 для текущего состояния ОУ; 


4) формирование сигнала оптимального управления (9). 

Для реализации этого алгоритма необходимо наличие адекватной 
математической модели нелинейного динамического процесса. 

Идентификация сложных нелинейных динамических ОУ традиционными 
способами требует больших расходов на экспериментальные исследования. Методы 
же нелинейной динамики позволяют с единых позиций определять 
(классифицировать) и исследовать режимы функционирования технологических 
процессов по отдельным временным реализациям, а также оценить структуру ОУ и 
синтезировать его модель. 

Нелинейный процесс (динамическая система) может быть описан [3] с 
помощью векторного уравнения потока или дискретного отображения Пуанкаре 
размерности 4 и содержащих вектор параметров порядка системы р. 


Динамические системы имеют в зависимости от значений параметров р 


четыре типа решения [3]: равновесие, когда после переходного процесса система 
достигает стационарного состояния, периодическое и квазипериодическое решения, 
а также хаос. Этим решениям соответствуют аттракторы системы в виде 
устойчивого равновесия, предельного цикла, квазипериодического и хаотического 
аттракторов. При изменении параметра р система теряет устойчивость своего 
состояния (режима функционирования) и переходит в другое состояние 
(бифуркация). 

Известно [3], что по одной временной реализации можно определить 
корреляционную энтропию, которая характеризует оценку глубины точного 
прогноза состояния порождающей системы и в каком режиме она находится, а также 
корреляционную размерность аттрактора (порядок системы). 

Энтропия Колмогорова К описывает динамическое поведение на аттракторе и 
пропорциональна скорости потери информации о состоянии динамической системы 
во времени. Для регулярного движения К -энтропия равна нулю, для систем с 
детерминированным хаосом — положительна и постоянна. Энтропия К бесконечна в 
случае поведения системы как белого шума, что говорит об отсутствии 
предсказуемости процесса. 

Для оценки К-энтропии по экспериментальным данным используется 
величина корреляционной энтропии [3]: 


т К 
К, = Ша Ши ш ея <К 
&—0 Ко В (=) 


(10) 
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где № - длина временной реализации; = -— размер ячейки фазового 
а. $ 
пространства; К, (=) = Шт == У -к-х|| ]= УР. — обобщенный 

т и. 

› тр 
корреляционный интеграл; Уе-|к- х| ] — число пар Г и 1], для которых 
1 Е 
5] 


, 
К 
— 2 ый — 
№ ступенчатая функция 
п=0 


расстояние |х-х| <2; | 


Хевисайда. 
Дополнительно К - энтропия позволяет определить среднее время, на которое 
можно предсказать состояние системы. При этом, точное предсказание возможно 


КТ, 
только на интервале времени Т„ таком, что &.е`” =1 [3], тогда оценка интервала 


предсказуемости по экспериментальным данным выполняется как: 
Т, [127 (11) 
Кпр К, Е тир? 

Корреляционная размерность аттрактора характеризует сложность аттрактора 
динамической системы, то есть характеризует минимальное количество отсчетов 
(глубину памяти) переменных, входящих в математическую модель системы. 

Расстояние между ближайшими точками аттрактора до и после бифуркаций 
находится в универсальном отношении. Самоподобие такого явления описывает 
фрактальная размерность Хаусдорфа Р. Ее оценка, также как и оценка К - энтропии 
может быть получена по экспериментальным данным. 

Численная оценка размерности Ш) осуществляется в соответствии с 


алгоритмом Грассбергера-Прокаччиа вычисления корреляционной размерности [3]: 


. 105 В(Е 
рр. (12) 
=—0 102 = 
оч б 
где №(=) — количество ячеек; А(=)= Пт ле. -|* —х И — корреляционный 
№М->® р 
1] 
интеграл, учитывающий вероятность того, что две точки на аттракторе лежат внутри 
ячейки =“ (вероятность того, что две точки разделены расстоянием, меньшим =). 
Для определения Д, строят зависимость 105К(=) от 105= и ищут на ней 
линейный участок, наклон которого и определяет искомое значение корреляционной 
размерности Д,, которое далее используется для определения размерности 


аттрактора. Размерность 4, начиная с которой корреляционная размерность О, 


перестает изменяться, является минимальной размерностью вложения аттрактора 
(наименьшей целой размерностью фазового пространства, которая содержит весь 
аттрактор). 

Вместе с тем, из теоремы про вложение [3] вытекает, что оценка размерности 
4 определяется через оценку размерности аттрактора ДО, как: 


ао. + (13) 
Таким образом, вычисление размерности ДР, позволяет реконструировать 


аттрактор, а также определить размерность 4 переменных математической модели 
ОУ (глубину памяти входных и выходных переменных модели). 
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Наличие у нелинейного динамического процесса различных режимов 
функционирования требует выполнения не только параметрической, но и 
структурной идентификации ОУ. 

Процесс структурно-параметрической идентификации включает определение 
структуры, оценку и оптимизацию параметров модели ОУ. Первые две операции 
решаются путем генерирования моделей-претендентов различной сложности 
(базисных функций) и настройки их параметров с последующей селекцией лучших 
из них по выбранным критериям (оптимальная структура). Операция определения 
оптимальных параметров решается известными методами оптимизации путем 
уточнения полученных ранее значений параметров по критериям регулярности на 
всей выборке исходных данных [4]. 

Задача идентификации ОУ формулируется таким образом: на основании 
экспериментального множества функций (временных рядов) возмущений, 
управлений и выходов в условиях помех определить структуру (обобщенную 
функцию Р) и вектор параметров а модели вида (1), которая достаточно точно (в 
смысле некоторого критерия) аппроксимирует ОУ относительно входных и 
выходных величин во всем функциональном пространстве. 

Актуальными здесь являются выбор базисных функций, в терминах которых 
осуществляется идентификация, выбор способа генерирования и селекции структур 
разной сложности (метода структурной оптимизации), а также выбор метода 
параметрической оптимизации и эффективных критериев селекции и оптимизации. 

Для аппроксимации функций перспективным является использование 
нейронных сетей (НС) и гибридных НС с нечеткой логикой, параметрами которых 
являются весы нейронов, а также коэффициенты их функций активации и 
принадлежности (для гибридных сетей). 

Базисная функция в виде НС прямого распространения (НСПР) со скрытым 
слоем представляется как [5]: 


Я+и] = УЕ {У [9-Е (У у [91-0], (14) 
0ЕР [ЕО тЕеО 
где Р — множество глубин памяти соответствующих входов; Е — 


активационная функция выходного слоя НС; О — множество входов нейронов; /— 
порядковый номер входа выходного слоя НС; у -— весовые коэффициенты 


выходного слоя; Ё, — активационная функция нейронов скрытого слоя; т- 


у 


порядковый номер входа НС; у, „ — весовые коэффициенты связи т-го входа и /-го 
нейрона; у, - вход НС. 

Параметрами настройки (обучения) этой НС являются {у,у„} са, . 

Входами НС (14) согласно (1) являются ис {у}, а ее 


структурными характеристиками — {Т,,Р,Е,,Н,’,,М „} с Е., где Т, — тип структуры, 


уз" [2152 


= О -размер скрытого слоя, М„ — метод параметрической оптимизации 


(функция обучения НС). К НСПР относятся персептроны, каскадные НС, вейвнеты и 
др. 

Базисная функция в виде гибридной НС с нечеткой логикой (Апй$) 
представляется как [5]: 


А +"]=»` У [9]. „А-0], (15) 


9еР тео 


142 © В.И. Корниенко, С.М. Мацюк 


ТУ5М 1561 - 5359. Штучний 1нтелект, 2015, № 1-2 


где 9] =И„ (9/9;  аык-в]= ПИ. -9]}; Ч =И(@,); 


г=[а,). 
Здесь И7 -— функция, обратная функции принадлежности промежуточного 


выхода т сети с параметрами аг› а, — значение промежуточного выхода; Ти- 


т 


произвольная {-норма моделирования логической операции «И»; [Г „- функция 
принадлежности нечеткого правила / входа т с параметрами ау. 
Параметрами настройки НС (15) являются {а,,а}са., а структурными 


характеристиками — {Т,Р,И„,Ё 


т Рр> х 


М} <Е,, где ’,сО - количество правил 
разложения по входам. 

На основе выполненных исследований предложена методика идентификации 
нелинейного динамического процесса, состоящая из этапов оценки состояния и 
характеристик ОУ и его структурно-параметрической идентификации. 

Оценка состояния и характеристик ОУ включает: 

а) процедуру определения режима функционирования ОУ по корреляционной 
энтропии КА, согласно выражению (10); 

6) процедуру оценки характеристик ОУ, содержащую: 

- вычисление корреляционного интервала предсказуемости (глубины 
прогноза) процесса Т„„ по выражению (11); 


- вычисление корреляционной размерности аттрактора ДО» согласно (12); 


- определение размерности вложения аттрактора 4 (размерности фазового 
пространства — глубины памяти) ОУ по графику зависимости 2, (4) и согласно (13); 


в) формирование оценки вектора состояния 2 согласно выражению (2). 

Структурно-параметрическая идентификация ОУ включает: 

а) формирование задачи идентификации, которая выполняется с выбором: 

- видов критериев структурной и параметрической идентификации; 

- базисных функций (НСПР согласно (14) и Апй$ согласно (15)); 

- метода структурной (глобальной) оптимизации (метода прямого случайного 
поиска (ПСП) и генетического алгоритма (ГА); 

- методов параметрической оптимизации (определяются выбором базисных 
функций); 

6) структурную идентификацию, которая осуществляется с помощью 
композиции методов глобальной и локальной оптимизации путем генерирования 
структур моделей-претендентов (базисных функций (14), (15) со своими 
структурными характеристиками) и селекцию лучших из них по критерию 
структурной оптимизации; 

в) параметрическую идентификацию, которая заключается в определении 
параметров модели оптимальной структуры путем ее обучения методом локальной 
параметрической оптимизации по критерию регулярности. 

Для оценки эффективности предложенных методики и алгоритмов было 
выполнено моделирование процедур идентификации и синтеза оптимального 
управления процессом крупнокускового дробления руд (ККД) на основе 
экспериментальных данных, полученных в условиях Ингулецкого горно- 
обогатительного комбината. 
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В результате расчетов для сигнала содержания класса +100 мм в дробленой 
руде Г. и определены корреляционные энтропия К,=0,405 и размерность 


р,=1,931. Для определения размерности фазового пространства вложения 


аттрактора (глубины памяти процесса) вычислялась ее оценка сверху по выражению 
(13), а оценка значения 4 снизу определялась по зависимости О,(4). Получено, что 


для решения задачи идентификации процесса ККД, порождающего сигнал Г, и, 


глубина точного прогноза составляет 4 такта, а глубина памяти - от 3 до 5 тактов. 

Установлено, что минимуму структурного критерия несмещенности [4] для 
процесса ККД отвечают базисные функции в виде каскадной НСПР (14), 
оптимизированных с помощью ГА. При этом количество нейронов в скрытом слое 
составляет 26, функция активации скрытого слоя — сигмоидальная, выходного слоя — 
линейная, алгоритм обучения НС -— метод Флетчера-Ривса (Е|еспег-Вееуе$). 

В качестве меры точности параметрической идентификации модели 
оптимальной структуры использовался критерий регулярности (относительной 
среднеквадратической ошибки), значение которого для модели процесса ККД 
составило 0,0357. 


Определенно, что для процесса ККД предлагаемое оптимальное управление с 
интеллектуальным прогнозированием обеспечивает в сравнении с существующими 
системами управления снижение в -2 раза ошибки управления и повышение 
производительности следующего в технологической линии процесса само- 
измельчения руд (за счет стабилизации содержания класса +100 мм в его входной 
руде) на 3-15 % (при коэффициентах вариации свойств руды от 0,1 до 0,8) со 
снижением удельных расходов электроэнергии. 
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КЕЗОМЕ 
В.И. Корниенко, С.М. Мацюк 
Оптимальное управление нелинейным динамическим процессом с 
интеллектуальным прогнозированием его состояния 

В данной статье разработан алгоритм синтеза оптимального управления 
нелинейным динамическим процессом по функционалу обобщенной работы. 
Использование синтеза по этому функционалу упрощает решение и приводит его к 
определению минимума на скользящем цикле управления с привлечением в реальном 
масштабе времени информации о состоянии объекта управления к новому циклу 
управления и его будущего состояния по прогнозирующей модели объекта. 

Разработана методика идентификации нелинейного динамического процесса, 
включающая этап оценки его состояния и характеристик, а также этап его структурно- 
параметрической идентификации. Она позволяет по величинам корреляционных 
энтропии и размерности процесса определять размерность фазового пространства 
(глубину памяти) и реконструировать модель режима функционирования процесса. 

Реализация оптимального управления с интеллектуальным прогнозированием 
позволяет решить задачу повышения эффективности управления нелинейным 
динамическим процессом в условиях возмущающей среды функционирования. 


У.1. Когшепко, 5.М. Макшк 
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В.Т. Корниенко, С.М. Мацюк 
Оптимальне керування неллнШним динам!чним процесом з штелектуальним 
прогнозуванням його стану 

У данй статт! розроблений алгоритм синтезу оптимального керування нелйним 
динам!чним процесом за функцюналом узагальнено! роботи. Використання синтезу за 
цим функщоналом спрощуе рипення 1 приводить його до вдшукання мимуму на 
ковзаючому цикл! керування 13 залученням в реальному масштаб! часу 1нформаци про 
стан об'екту керування до нового циклу керування 1 його майбутнього стану по модел! 
об’екту, що прогнозуеться. 

Розроблена методика 1дентификаци нелийного динам!чного процесу, що включае 
етап ощнки його стану 1 характеристик, а також етап Його структурно-параметричнот 
1дентифткаци. Вона дозволяе за величинами кореляшиних ентропй та розмрност! процесу 
визначати розмрнеть фазового простору (глибину пам’ят!) 1 реконструювати модель 
режиму функцюнування процесу. 

Реалзащя оптимального керування з 1нтелектуальним прогнозуванням дозволяе 
вирипити задачу шдвищення ефективност! керування нелинйним динам!чним процесом в 
умовах збурюючого середовища функщонування. 
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